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Abstract 

Wc  study  fundamental  issues  in  electromagnetic  scattering  theory, 
with  an  emphasis  on  pole  behaviors  of  a  lossless  sphere  arising  from 
the  singularity  expansion  method  (SEM).  Wc  use  Mie  Theory  to  solve 
the  acoustic  and  electromagnetic  scattering  problems  for  spheres  with 
lossless  boundary  conditions  and  an  incident  plane  wave.  Wc  show 
that  for  certain  lossless  sheet  impedance  boundary  conditions  there 
exist  second  order  poles  for  both  cases.  Our  general  procedure  to 
directly  construct  lossless  sheet  impedance  boundary  conditions  which 
will  produce  high  order  poles  is  discussed  as  well  as  the  difficulties 
to  which  it  leads.  In  the  electromagnetic  scattering  case,  Foster’s 
Theorem  is  imposed  on  the  impedance  condition  to  ensure  that  a 
lossless  scattering  problem  is  obtained.  Wc  also  study  the  validity  of 
the  forward-scattering  theorem  associated  with  SEM. 
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1  Introduction 


The  singularity  expansion  method  (SEM)  [2]  was  introduced  in  1971  as  a 
way  to  represent  the  solution  of  electromagnetic  interaction  or  scattering 
problems  in  terms  of  the  singularities  in  the  complex-frequency  (s  of  two- 
sided-Laplace-transforrn)  plane.  Particularly  for  the  pole  terms  associated 
with  a  scatterer  (natural  frequencies),  their  factored  form  separates  the  de¬ 
pendencies  on  various  parameters  of  the  incident  field,  observer  location, 
and  scatterer  characteristics,  with  an  equally  simple  form  in  both  frequency 
(poles)  and  time  (damped  sinusoids)  domains. 

The  forward-scattering  theorem  is  a  classical  result  in  electromagnetic 
theory.  A  recent  paper  [4]  has  generalized  the  forward-scattering  theorem 
with  particular  application  to  lossless  bodies.  From  the  forward-scattering 
theorem  we  have  relations  between  the  absorption  and  scattering  cross  sec¬ 
tions,  and  the  forward  scattering.  The  scattered  fields  are  represented  by  a 
scattering  dyadic  times  the  incident  plane  wave.  This  allows  one  to  refor¬ 
mulate  the  results  in  terms  of  the  scattering  dyadic,  exhibiting  some  general 
characteristics  of  this  dyadic.  It  is  extended  (for  lossless  scatterers)  by  ana¬ 
lytic  continuation  away  from  the  jus  axis  out  into  the  complex  s  plane  and 
applied  to  poles  in  the  singularity  expansion  method.  In  particular  this  gave 
new  insight  into  the  scattering  natural  frequencies  and  modes,  also  implying 
new  ways  to  calculate  them  from  the  scattering  operator  in  the  right  half  of 
the  complex  s— plane.  This  gives  new  insight  into  the  properties  of  the  poles 
(natural  frequencies,  sa)  in  the  left-half  plane  and  the  associated  natural 
scattering  modes. 

Although  in  practice,  we  only  encounter  the  first  order  scattering  poles,  an 
interesting  question  concerning  the  SEM  and  the  forward-scattering  theorem 
concerns  the  existence  of  higher  order  scattering  poles.  Carl  Baum  showed 
that  2nd  order  poles  can  be  constructed  for  a  transmission  line  problem  [6]. 
Since  the  transmission  line  problem  is  finitely  dimensioned,  we  can  actually 
use  the  scattering  matrix  to  find  the  poles  (i.e.  the  eigenvalues).  However,  in 
general  the  problem  is  infinitely  dimensioned.  Thus,  we  consider  a  classical 
model  problem,  scattering  from  a  sphere  with  an  incident  plane  wave.  We 
compute  the  exact  solution  using  Mie  Theory.  In  [2],  Carl  Baum  showed  that 
for  a  perfectly  conductor  sphere,  there  only  exist  first  order  poles.  Saucer 
12]  also  proved  some  similar  results  for  a  general  shape  scatterer. 

We  show  that  for  the  acoustic  scattering  problem  high  order  poles  can 
be  constructed  for  certain  impedance  boundary  conditions,  while  for  hard 
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and  soft  spheres  there  only  exist  first  order  scattering  poles.  The  general 
procedure  to  construct  arbitrary  order  poles  is  discussed  as  well  the  nec¬ 
essary  condition  for  the  conservation  of  energy  to  the  scattering  problem. 
For  electromagnetic  scattering  from  a  sheet  impedance  loading  sphere,  we 
show  that  there  exist  2nd  order  poles.  Foster’s  Theorem  is  imposed  on  the 
impedance  condition  to  ensure  the  scatterer  is  lossless.  The  existence  of  a 
2nd  order  pole  for  a  transmission  line  scattering  problem  is  included  in  sec¬ 
tion  4.  Some  analytical  and  numerical  treatment  of  the  forward-scattering 
theorem  is  presented  in  section  5. 


2  Scattering  from  a  lossless  acoustic  sphere 

2.1  Introduction 

In  this  section,  we  are  considering  the  following  problem.  An  incident  plane 
sound  wave  is  propagating  in  some  direction  with  the  scatterer  being  a  sphere. 
The  scattered  solution  can  be  written  explicitly  using  spherical  harmonics. 
However,  the  thing  we  are  really  interested  in  is  to  explicitly  solve  the  prob¬ 
lem  using  SEM  and  to  study  the  scattering  pole  behavior  of  the  solution. 
By  putting  different  boundary  conditions  on  the  scatterer,  we  are  able  to 
construct  not  only  simple  poles  but  also  2nd  or  even  higher  order  poles. 
Meanwhile,  we  always  keep  the  impedance  function  satisfying  the  Foster's 
Theorem,  so  that  it  remains  a  lossless  system. 

2.2  Formulation  of  the  acoustic  scattering  problem 

Considering  the  linear  acoustic  equation  [11] 


^  +  PoV-,-0 

(1) 

dv  _ 

A,s+Vp=0 

(2) 

p  =  c2p 


c2  = 


0 


where  p  is  for  acoustic  pressure,  p  is  for  density,  v  is  for  fluid  velocity,  c 
is  the  speed  of  sound.  We  want  to  solve  for  the  scattering  solution  to  the 
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Figure  1: 
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above  equation  with  an  incident  plane  wave  and  some  different  but  lossless 
boundary  conditions  on  the  sphere.  We  will  use  the  expansion  in  terms  of 
spherical  harmonics  to  solve  the  equation  in  order  to  locate  the  poles.  In 
fact,  this  is  just  the  SEM  m  the  acoustic  case.  Let’s  first  derive  the  wave 
equation  from  the  system  above,  so  that  later  on  we  will  treat  the  problem 
mathematically,  temporarily  disregarding  its  physical  meaning.  Plugging 
p  =  c2p  into  (1)  we  get 


±+n>c2V-v  =  tt 

(3) 

8v 

p*dt +  Vp  - 0 

(4) 

if  we  ^(3)  and  plug  in  (4),  we  get 
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similarly  for  v,  we  ^  (4),  and  plug  in  (3)  we  get 


S  +  —  V(-p„c2V  •  v)  =  0 

OT  p0 


d2t> 


-  c2V(V  •  u)  =  0 


assume  v  is  irrotational  i.e.  V  x  v  =  0,  then  V(V  •  v)  =  V2u.  We  get 


d2v 

dt? 


-  c2V2r  =  0 


Thus  if  u  is  either  p  or  v, 


(6) 


d2u 

w 


-  c2V2u  =  0 


Taking  the  Laplace  transformation  we  get 


(V2  —  7  2)u  —  0,  where  7  =  -  (7) 

c 

Assume  the  scatterer  is  a  ball  with  radius  1  and  impedance  boundary 
condition  of  the  following  form 

Oil 

—  +  a(s)u  =  0  ,  at  r  —  1  (8) 

Oil 

where  n  is  the  outward  normal. 

Assume  the  incident  plane  wave  is  in  the  2-direction  i.e. 

{t(inc)  =  e-7(0,O,l)  (x,2/,z)  _  e~7Z  ^ 

We  want  to  solve  scattered  solution  explicitly  according  to  (7)(8)(9) 
and  study  the  pole  behavior  of  the  scattered  solution  u^scK 


2.3  Hard  and  soft  spherical  scatterer 

First  of  all,  we  want  to  relate  the  mathematical  impedance  function  a( 7) 
to  the  actual  acoustic  impedance  Za(s).  Units  of  acoustic  impedance  are 
Pa  ■  s/vi  or  kg/(m2  ■  s)  By  definition 


Za(s) 


T>(r,  fi) 
v(r , s)  •  nlH 


w  here  p,  v  are  the  Laplace  transforms  of  p  and  v,  and  nul  is  the  inward 
normal.  Take  the  Laplace  transform  of  (2),  we  get  sv+  ^Vp  =  0.  And  then 
take  the  inner  product  with  outward  normal  direction  n  ,  we  get 


1  dP  „ 

sv  •  n  H - —  =  0 

Po  an 

dp 

=>  Q-  =  -Pqsv  •  n 
an 

Then  acoustic  impedance  can  be  written  as 


(10) 


Z(s) 


P 


v  •  n. 


P 

^pos 

dn 


Iii  our  formulation,  if  we  plug  (10)  into  (8),  assuming  we  put  the  impedance 
oil  pressure(i.e.  u  —  p  here),  we  derive 


— —  =  — ru  =  z(s 

v  •  nin  o(.s) 

which  relates  the  mathematical  impedance  condition  to  the  real  acoustic 
impedance  boundary  condition. 

The  infinite  specific- acoustic-impedance  limit  \Z\  — *  oc  corresponds 
to  a  hard(rigid)  surface,  the  limit  \Z\  — >  0  corresponds  to  a  soft  (pressure- 
release)  surface.  Thus,  for  a  hard  sphere  the  mathematical  impedance  bound¬ 
ary  conditions  become  =  0  and  u  =  0  for  the  soft  sphere. 

For  detailed  computation  of  the  following  please  see  section  5. 


{l(mc)  =  e  yZ  =  e-7rcos  0  =  ^(2«  +  1)  (■ -  1 )"  ln  (7r)  P„(cOS  6) 

n— 0 
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do  rn=n 

U {sc)  =  E  E  «nm(7)Aw(7r)>;r(0^) 

n=0  m=— n 

DO 

=  E  an(7)^n(7r)^rt(CQS  0) 

n=0 


Where  kn(s)  is  the  modified  Bessel’s  function,  an( 7)  is  a  coefficient  to  be 
determined. 

Applying  the  impedance  boundary  condition  for  hard  and  soft  spheres 
respectively,  at  r  =  1 


dii 


(sc) 
ha  rd 


On 


dii 


(inc) 

hard 


On 


~  (sc)  ^  (inc) 

Usoft  =  ~Usoft 

n  is  the  outward  normal,  i.e.  r  in  our  case(spliere). 
We  derive  the  scattered  solutions  as  follows 


d«c) 

lhard 


E 


n=0 


(-ir1(2n+l)K(7)] 

7^(7) 


kniri')  Pnicos  0) 


d*c) 

lsoft 


E 

H=0 


( —  l)Tt+1  (2/?.  +  l)fon(7)] 

M  7) 


kn(^r)Pn(cose) 


Since  the  modified  spherical  Bessel  functions  kn  and  their  derivatives  only 
have  simple  zeros  [2],  for  both  the  hard  and  soft  sphere  scatterer,  the  solution 
only  has  simple  poles. 


2.4  Lossless  impedance  loading  of  a  sphere 

Consider  the  total  energy  in  the  usual  mathematical  way 

Rd  =  J  j  J 2  (ui  +  l^Ji|2) 
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dEo 

nr 


ututt  +  v/  •  VU( 


ut(utt  -  V2u 


dD 


and  the  net  energy  flnx  is 


dr 


(11) 


dD 

where  in  our  case  D  is  a  sphere.  Plugging  in  our  impedance  boundary  con¬ 
dition,  the  net  energy  flux  becomes 


~J  Jsa(s)\u(y,s)\ 2 

dD 

For  a  lossless  case,  the  total  energy  change  should  be  zero,  which  implies 
J1^(11)=0=>(11)  should  be  odd  =>  a(s)  should  be  even,  Z(s)  odd.  Math¬ 
ematically,  this  condition  is  necessary  to  guarantee  that  energy  is  conserved. 
However,  in  order  to  extend  our  results  later  on,  we  need  more  constraints 
on  Z(s)  or  equivalently  on  a(s).  Namely,  we  want  Z(s )  to  satisfy  the  Foster 
reactance  theorem  which  we  will  discuss  a  little  bit  later. 


2.5  High  order  poles 
2.5.1  2nd  order  poles 


The  general  expansion  (see  section  5  for  more  details  about  the  calculations) 
of  the  scattered  solution  can  for  our  impedance  condition  can  been  written 

as 


„(“> = £ 
n= 0 


(— l)n+1(2n  +  l)[a(s)iw(a)  + 
a(s)kn(s)  +  sk'n(s) 


kn(sr)Pn(cosO) 


Assume  c  =  1  here,  so  7  =  s 
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The  far  field  pattern  of  the  solution  is 


u<“>  = 


sr  ~ 

— £ 

n—0 


(— l)n+L(2n  +  l)[g(s)in(s)  +  ^,,(5)] 
s(a(s)kn(s)  +  sk'n(s)) 


Pn(cos9) 


Ill  order  to  construct  a  second  order  pole  we  need  both  the  denominator 
and  the  derivative  of  the  denominator  for  the  scattered  solution  to  be  zero 
at  some  specific  s . 

That  is,  the  denominator^) 

sa(s)kn(s)  +  s2^-kn  ( s )  =  0  (12) 

as 

and  the  derivative  of  the  deriominator=0 


Oc  (s)  kn  (s)+s  ^  — a  (5)^  kn  (s)-t-sa  (s)  —  kn  (s)+2  s — kn  {s)-\-s2-j-^k :n  ( s )  -  0 

(13) 

Note  that  we  do  not  want  to  solve  the  above  system  of  ODEs,  because 
we  only  need  the  equations  to  hold  at  one  specific  s  for  some  pre-chosen 
impedance  boundary  condition  a(s) 

By  using  Bessel’s  equation  itself  we  can  replace 

s2~kn  (s)  =  2  s~j^kn  (s)  -  (,S2  +  n  ( 77.  +  1))  kn  ( s ) 

Solve  (12)  for  ^ kn(s )  and  plug  into  (13),  we  get 

k„  ( s )  (s  a  (,s)^  +  a(s)  ~  rv  Gs)~  +  s*"  +  n  +  n)  =  0 

for  kn(s)  ^  0  we  want 

sa'  +  s2  +  a  +  n(n  +  1)  —  a2  =  0 

Choose  n  =  l,k\(s)  =  f 
Solve  (12)  for  a(s),  we  get 

,  .  s2  +  2  s  +  2 
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Solve  (13)  for  a'(s),  we  get 


a'(s) 


s  (s  +  2) 
s2  +  2  s  +  1 


Note:  We  can  not  differentiate  a(s)  here,  since  it’s  just  the  value  evaluated 
at  s,  not  a  function. 

The  solution  we  derived  above  simply  means  that,  if  we  want  to  construct 
a  2nd  order  pole  at  s,  the  impedance  function  a(s)  should  be  chosen 


1.  a (s)  Even  function 

2.  Value  of  a(s)  at  s  equals 

3.  Value  of  a'(s)  at  s  equals 


Let’s  build  a  concrete  example  to  verify  the  above  results. 
Choose  s  —  —2,  then  a(s)  =  —2,  af(s)  =  0 
Assume  a(s)  has  the  following  form 


a(s) 


Ci  -I  c2.s2 
1  +  C3S2 


after  some  calculation  we  get  a(s)  =  —2 

The  denominator  of  the  scattered  solution  ii^  for  n  =  1  is 

—2sk\(s)  +  s2k[(s)  =  ^  ^  e-* 

2  s 

Thus,  we  have  a  2nd  order  pole  at  s  =  —2. 


2.5.2  3rd  order  poles 

The  process  is  similar  to  the  2nd  order  pole  case.  We  want  denominator, 
denominator denominator”  =0  at  some  s. 

sa(s)  kn(s)  +  s2-j-kn(s)  =  0  (14) 

as 


<*  (■s)  kn  (s)+.s  (  —a  («s)  )  kn  (,s)+«sa  (s)  —  kn  (.s)+2,s  —  kn  (•"i)+fi2^^k„  (a)  =  0 


ds 


ds 


ds 2 


(15) 
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2  (s)^  kn  (s)  +  2a  (s)  -^kn  (s)  +  s  («)  j  k,  (s)  +  2s  (s)^  ~ku  (s) 

+  sft(s)  -^A :„  (s)  +  2  ~kn  (s)  +  4s-^kn  (s)  +  s2  -^kn  (s)  =  0  (16) 


B\  computing  the  derivative  of  Bessel’s  equation  we  get 


2  -kn  (s)  =  -4s^kn  {s)-2  -^-kn  (s)+2  skn  (s)+(s2  +  n(n  +  1))  -^ku  (,s) 


ds3 


together  with 


s 


2 


kn  (s)  -  (s2  +  n  (n  +  1))  kn  (s) 


and 


d  /„\  _  <*(S)M«) 

Kn[S)  — 

as  s 


equation  (15),  (16)  reduce  to 


d 


k„  (s)  (s)  +  s~nQ  00  -  (Q  (s))2  +  s2  +  7)2  +  71 

k„  (s)  (2  £q  (,)  +  s^a  (s)  -  2  (J-sa  (s))  a(s)  +  2s 
solve  (14),  (17),  (18)  for  a,  a', a"  respectively,  we  get 


=  0 
=  0 


a(s)  =  — 


(  ■>  ) 

kn  (s) 


(17) 

(18) 


o'(s)  =  - 


~s  (s))  kn  (s)  ~  S2  (^A’n  (s))  +  6~2  (kn  (s))2  +  1 12  (fc„  (s))2  +  II  (A~„  (.s))~ 

~T{kj7)f~ 
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a"(s)  =  2{s(£kn(s)')(kn(s))2-2kn(s)s2^kn(s^  +  n2  (kn  (s))3 

+  n  (kn  ('S))J  -  (s)^  +  S3  (lS)^  (kn  («))2 

+  s  (j^kn  (s)^  n2  (kn  (s))2  +  s  (j^ku  (s)^  n  (kn  (s))2}/s2  (k„  (s))3 

Choose  n  =  2,  /q>(s)  =  \  — — ^3——^,  we  get 


q(s) 

«'(a) 


a"(s) 


s3  +  4s2  +  9s  +  9 
s2  +  3s  +  3 
s(s3  +  6  s2  +  12s  +  6) 
(s2  +  3s  +  3)2 
6(3  +  s3  +  6s2  +  9s) 
(s2  +  3s  +  3)3 


Again,  the  solution  we  derived  above  simply  means  that,  if  we  want  to 
construct  a  third  order  pole  at  s,  the  impedance  function  o(s)  should  be 
chosen 


1.  a(s)  Even  function 

2.  Value  of  o-(s)  at  s  equals 

3.  Value  of  ot'(s)  at  s  equals 

4.  Value  of  a"(s)  at  s  equals  r^'{^/+gs6^3|39s^ 

Let  test  our  results  again  with  s  =  —4.  Thus,  we  require  a(— 4)  =  —  y, 

“'(-4)  =  i>  «"(-4)  = 

Assume  a(s)  still  has  the  following  form 


q(s)  = 


C\  +  C2S2 
1  +  c3s2 
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we  get 


The  denominator  of  the  scattered  solution  for  n  =  2  is 


7 re  s  T  3 s  +  12)  (s  +  4) 
2  s2(54  +  s2) 


Thus,  we  have  a  3rd  order  pole  at  s  =  —4. 

2.5.3  Aribitrary  order  poles 

In  principle,  we  can  follow  this  procedure  to  get  poles  of  any  high  order  at 
any  specific  location.  However,  the  computation  will  become  messier  and 
messier.  It  will  become  more  clear  to  construct  high  order  poles  if  we  try  to 
view  this  process  using  Taylor  expansion  around  the  pole  sp. 

First,  let/s  introduce  some  very  nice  properties  of  the  modified  Bessel 
function  [9]. 


where,  for  /  ^  0, 


/-i 


(2/  -  A’)! 


We  also  have  the  following  continued  fraction  representation  for  S\\ 


Thus,  it  is  possible  for  us  to  rewrite 
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*,(.) = - 

sf  (»)  s 

where,  /c;ll)(s),  /c;2)(s)  are  just  polynomials  in  s. 

Suppose  impedance  function  a(s)  has  the  following  form 


o(s)  =  -cq- 


(s~  +  ao)(s2  +  02)  •  •  •  (s2  +  at) 


-Co- 


ai(s) 


(s2  +  ai)(s2  +  a3)---(s2  +  a2„-i)  a2(s) 

where,  Co  >  0,  0  <  a0  <  a\  <  a2  <  •  •  •  <  a 271-2  <  ^2n-i  <  ^2?i  <  oc, 
A;  —  2n  —  2or2n,  and  ai(s),  a2(s)  are  of  course  polynomials. 

For  example,  for  2nd  order  pole  we  need  2  free  parameters,  thus  a(s)  can 
been  chosen 

cv(s)  =  -Co(s2  +  a0) 

For  3rd  order  pole  we  need  3  free  parameters,  thus  o(s)  can  been  chosen 


a(s)  =  -c0 


(s2  +  a0) 


( s 2  +  ai) 

For  4th  order  pole  we  need  4  free  parameters,  thus  a(s)  can  been  chosen 

(s2  +  a0)(s2  +  o2) 


a 


(s)  =  -CQ- 


(s2  +  fll) 

The  denominator  to  the  scattered  solution  can  now  be  written  as 
kn\s)  e~s  (  ot\ (s) 


(co 


a2(s) 


+  s  +  1  +  S„(s 


(19) 


k(n\s)  S 

If  we  want  to  construct  a  j— th  order  pole  at  sp,  we  just  need  to  choose  the 
correct  power  of  a\ (s),  <a2(s)  and  let  n  =  j  —  1,  and  rewrite  the  numerator 
using  Taylor  expansion  around  sp. 


j- 1 


Hi\s)(^c0  <*l{s)Qn(s)  +  a2  ($)f  u(*s))  “  (Ms  ~  spY  +  -  spY)  (20) 

i=0 


where  s  +  1  +  Sn(s)  = 

Wn\S) 
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Choose  Co,  a0,  ai,  . . .  aj_2  so  that  fa  =  0  for  i  =  0, . . .  j  —  1.  Note  that 
we  always  want  Co  and  aj  satisfy  our  assumptions. 

Thus  we  construct  a  lossless  impedance  function  a($),  which  will  produce 
a  j— tli  order  pole  at  sp.  One  might  ask:  Can  we  always  get  a  solution  which 
satisfies  all  of  our  assumptions  for  arbitrary  sp  <  0?  At  least,  in  our  case, 
the  answer  is  NO!  There  will  be  some  restriction  on  sp. 

Let’s  go  over  the  3rd  pole  again  using  the  Taylor  method  to  get  a  better 
picture  of  the  procedure  and  the  restriction. 

We  are  constructing  3rd  order  pole  at  s  =  sp,  so  j  =  3,  n  =  2, 


a(s) 

h{s)  = 


(s2  +  a) 

~  _c(s2  +  6) 

e~s  ( s 2  +  3s  +  3) 


s  T  1  +  ^(s)  — 


Pn(s)  s3  +  4  s2  +  9  s  T  9 


Qn(s)  s  +3s 
According  to  (20),  the  zeros  are  contained  in 


(s2  +  3 s  4-  3)^c(s2  +  a)(s2  4-  3  s  +  3)  +  (s2  +  b)(s 3  4-  4s2  4-  9  s  4-  9 (21) 
expand  (21)  around  (s  —  sp ),  we  get  the  coefficients 
0o  = 

sp‘  T  18 csp3  b  7 bspA  4“  bsp 3  +  csp ^  4-  6 casp3  4-  9c«Sp2  +  24 sp^  4-  7 sp()  b 
15  cdSp 2  4~  casp A  +  48  spA  +  6  csp  3  4”  54  sp3  4"  24  bsp 3  4~  15  csp ^  +  54  bsp  +  18  casp  + 
27  sp2  +  27b  +  48  bsp2  4-  9  ca 

0i  = 

54  csp2  b  7  Sj/  4"  30  casp  +  4  ccisp3  +  6  csp 3  -b  162  sp2  4*  18  ccisp2  b  120  sp  *  + 
30  csp  ^  T  18  cci  T  72  bs p2  +  28  bsp3  4"  18  csp  +  42  sp  y  T  54  sp  +  96  bsp  4"  192  sp3  b 
5  bspA  4~  54  b  +  60  csp3 

02  = 

72  bsp  4-  6  ccisp~  +  162  sp  +  54  csp  +  15  cspA  +  18  casp  T  21  sp^  +  105  sp  *  + 
27  +15  ca  +  60  cSp3  +  48  6  +  240  sp3  +  9  c  +  10  bsp 3  +  288  sp2  +  90  csp2  +  42  bsp2 


15 


Solve  po  =  /3\  =  (32  =  0  for  a,  b.  c  in  terms  of  sv  we  get 


+  12  Sp 1  +  81  Sp^  +  315  Sp ^  +  648  Sp^  +  648  Sp  +  216 
3  6’/  +  28  sp3  +  99  sp2  +  153  sp  +  96 

b  g  ”b  ^  Sp ^  +  35  Sp^  +  72  +  /  2  Sp  +  24 

sp3  +  9  Sp 2  +  27  Sp  +  24 


3  Sp4  +  28  3p^  +  99  Sp' ^  +  153  Sp  +  96 
sp3  +  9  5p2  +  27  sp  +  24 


If  we  choose  sp  =  —4,  then  a  =  yy,  ft  =  54,  c  =  11,  clearly  c  >  0,  and 


0  <  a  <  b  <  oc 


a(s) 


-H(«2  +  fj) 

s2  +  54 


which  is  the  same  as  we  showed  before.  In  fact,  if  we  choose  sp  —  —3, 
then  a,  6,  c  will  not  satisfy  our  assumptions.  If  we  plot  a,  b ,  c  in  terms  of  syj, 
we  will  see  that  in  order  to  satisfy  all  the  assumptions,  sp  can  only  be  chosen 
approximately  sp  <  —3.2.  At  sp  =  —3,  a  will  be  negative.  For  the  more 
general  case,  we  can  also  use  this  graph  to  determine  the  range  of  sp.  For 
example,  in  this  case,  the  range  of  sp  will  be  the  interval  where  the  graph  of 
b  is  above  a  and  all  of  the  graphs  a,  b)  c  are  above  the  x-axis. 


2.6  Interpretation  of  the  Results 

Our  prototype  here  for  scattering  is  a  constant-frequency  plane  wave  pro¬ 
ceeding  in  direction  ez.  The  overall  acoustic  pressure  is  written 

p(s,  r,  9)  =  p(mc)(s,  r,  6 )  +  p^c)(s,  r,  6) 

where  p^inc\s}  r,  9)  is  the  incident  plane  wave  and  p^sc\s,  r,  9)  is  the  scat¬ 
tered  solution( wave’s  complex  amplitude).  The  function  p(sc\s,  r,  9)  satisfies 
the  Helmholtz  equation  and  the  Sommerfeld  radiation  condition. 

For  a  hard  surface  scatterer  we  require 

Vp  =  0  =>  Vp(5c)  •  n  =  - Vp(inc)  •  n 

For  a  soft  surface  scatterer  we  require 
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Figure  2:  plots  of  a.b,c 
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p  =  0 


p(s°) 


n  =  -p(inc) 


n 


where  n  is  a  unit  normal  vector  pointing  into  fiuid. 

If  we  assume  the  scatterer  is  a  sphere  centered  at  the  origin  with  radius 
R.  We  have 


-(sc) 

Phard 


OC 

-E 

n=0 


(-l)"+1(2n  +  l)K(7fl)] 
7^(7  «) 


kn(^r)Pn(cos0) 


-(sc) 

Psoft  =  2 L 


(-l)"+1(2»+1)[7»n(7)] 

kn(l  R) 


k„(^r) Pn(cos  0) 


The  poles  of  p\l^ld  sojt  are  all  simple  as  discussed  above.  Since  vl'sc'1  have 
exactly  the  same  zeros  as  p*scMoes,  we  will  just  talk  about  7/  5C1  here. 

If  the  sphere  scatterer  has  a  specific-acoustic-impedance  (neither  soft  nor 
hard)  condition  (equivalent  to  rv(.s)  =  —2) 


then 


Z(s)  = 


po£ 

2 


rPL' 


=  E 


Tf  —0 


(-l)w+1(2»-  +  l)[a(7K)*n(7ft)  +  7^n(7/Q] 

<*W)knW)  +  yRk'niiR) 


kn  (77 )  Pn  (COS  0) 


Thus  p(*c)  will  have  a  second  order  pole  at  7/i  =  —2  or  s  = 

If  the  sphere  scatterer  has  a  specifie-acoustie-iinpedanee  (neither  soft  nor 


hard)  condition  (equivalent  to  a(s)  = 


zHff!±ix2 

*2+54 


) 


_  Posjs2'  +  54) 

(s)  n(»2+?t) 

then  p^c)  will  have  a  third  order  pole  at  7 R—  -4  or  s  =  — 

I11  12],  Saucer  showed  that  for  acoustic  scattering  SEM  poles  are  simple, 
but  there  is  no  contradiction.  I11  that  reference,  the  impedance  condition 
is  not  considered  in  the  proof  and  Sancer  used  a  hard  acoustic  scatterer  to 
derive  the  simple  pole  behavior  for  arbitrary  shape  of  scatterer,  with  which 
our  results  agree. 
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3  Scattering  from  a  lossless  electromagnetic 
sphere 

3.1  Introduction 

In  this  section,  we  are  considering  the  problem  of  a  plane  wave  incident  on  a 
sphere  (with  perfectly  conducting  surface  and  lossless  sheet  impedance  load¬ 
ing  respectively)  as  illustrated  in  figure  3.  An  E  wave  has  been  chosen  as  an 
incident  electromagnetic  plane  wave  propagating  in  l!  direction.  The  scat¬ 
tered  solution  as  well  as  the  surface  current  density  can  be  written  explicitly 
using  vector  spherical  harmonics.  In  the  case  1,  a  perfectly  conducting  sur¬ 
face,  we  will  briefly  summarize  the  work  done  by  Carl  Baum,  which  shows 
that  there  exist  only  first  order  scattering  poles.  In  the  case  2,  a  lossless 
sheet  impedance  loading  sphere,  there  exist  2nd  order  scattering  poles  for 
some  mathematically  chosen  boundary  conditions.  Foster’s  Theorem  is  en¬ 
forced  on  the  impedance  function  Zs(s )  to  guarantee  it  is  a  realizable  physical 
boundary  condition. 

3.2  Formulation  of  the  electromagnetic  scattering  prob¬ 
lem 

Define  a  set  of  orthogonal  (right-handed)  unit  vectors  by 

U  =  sin(<91)cos(01)lx  +  siii(01)sin(01)fy  +  cos(9l)lz 

12  =  -  cos($! )  cos(0! )  lx  -  cos(^i)sin(01)ly  +  sin(0i)lz 

13  =  sin(0i)lx  —  cos(0i)ly 

As  shown  in  figure  4,  l\  is  the  direction  of  propagation  and  12  and  13  are  mu¬ 
tually  orthogonal  unit  vectors,  each  orthogonal  to  l\  to  indicate  the  polariza¬ 
tion  of  the  electromagnetic  fields  in  the  incident  plane  wave.  For  convenience 
12  is  chosen  in  a  plane  parallel  to  li  and  the  z  axis  (E  or  TM  polarization  if 
the  electric  field  is  parallel  to  12)  while  I3  is  parallel  to  the  x)  y  plane  (H  or 
TE  polarization  if  the  electric  field  is  parallel  to  I3).  In  free  space,  electro¬ 
magnetic  plane  waves  have  both  electric  and  magnetic  fields  orthogonal  to 
li.  Thus  only  12  and  I3  are  concerned.  This  removes  the  L  functions  (details 
are  shown  later)  in  the  expansion  (plane  waves  have  zero-divergence  fields). 
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Figure  3: 


sheet  impedance  Zt  ( s ) 


Figure  4: 


21 


We  can  use  the  relations  between  Cartesian  and  spherical  coordinates 

x  =  r  sin(0)  cos(0) 
y  —  r  sin(0)  sin(0) 

2  =  V  COS(0) 


lx  =  sin(0)  cos(0)l,  +  cos(0)  cos(0)4  —  sin(0)4 
I y  =  sin(0)  sin(0)lr  +  cos(0)  sin(0)4  +  cos(0)4 
4  =  cos(0)4  —  sin(0)4 

to  express  the  incident- wave  unit  vectors  in  terms  of  and  (0. 0)  as 

li  =  [cos(^i)  cos(0)  +  sin(0i)  sin(0)  cos(<^>  —  <4)]4 
+  [—  cos(4)  sin(0)  +  sin(4)  cos(0)  cos(<p  —  <4)]4 
+  [-  sin(0x)  sin(<p  -  <£i)]4 


4  =  [sin(0i)  cos(0)  —  cos(0i)  sin(0)  cos(<^)  —  0i)]lr 
[sin(^i)  sin(0)  +  cos(4)  cos(0)  cos(<?!>  -  <4)]4 
+  [cos(0i)siii(</>  -  <4)]4 


4  =  —  sin(0)  sin(0  —  <4)4 

—  cos(0)  COS($  —  <4)4 

-  cos(<£-<4)]4 


Having  the  direction  of  incidence  and  two  polarizations  expressed  in  spher¬ 
ical  coordinates  we  can  go  on  to  express  the  response  to  some  delta  plane 
wave  functions.  For  an  incident  delta  function  plane  wave  we  need  spherical 
harmonics  and  vector  wave  function  in  which  to  express  the  expansion  in 
spherical  coordinates.  In  spherical  coordinates  we  have  the  common  differ¬ 
ential  operators  as 


VF  - 
V  ■  F  — 
V  x  F  = 


d  _  -  1  d  „  -  1  d 

1  r22~F  +  lg-  — qF  +  4  •  7Z\  a  1 1 

dr  rod  rs\n{d)d(p 


j__a,2  1  d 

v2  dr  '  r  sin(0)  86 


(sin  {9)Fe)  + 


1  d 

r  sin(0)  d(p 


d 


S'  sin (6)  d6 


(sin (0)F0)  - 


d 


7'sin(0)  d(p 


VSF 
Vs  ■  F 
V,  x  F 


*'WF  +  i* 


1  0 

sin(0)  d(j) 


F 


1  d 

sin (0)  dO 

-  .  1 
lrh 


1 


lsin(0)  dO 


(Sm(«)F,)  +  ,if  , 

9  (sin(6)Ft)-~~l-,\  +  U 


sin  (6)  d<p 


1  d 
sin(0)  d<j) 


Fr 


Spherical  Harmonics 

The  scalar  spherical  harmonics  can  be  written  as 


Yn,m'(e,<t>)  =  P}r)(cos(0)) 


cos(m(f)) 

sin(m<f>) 


where  Pfin\x)  is  the  Legendre  function  defined  as 


rfm 

Pn(x)  =  -  x2r^—Pn(x),  Pn(x)  =  f#{x)  = 

Vector  spherical  harmonics  are  defined  as  follows 


1  dn 
2 nu\  dxn 


(x2  -  1) 


4*.m,p(0,  <j>)  =  y^niA®!  0)  lr 

0)  —  0)  lr  X 

^n.m,p(^)  0)  X  /'n,m,p(^i  0)  —  lr  X  Qn,rn,p 


They  also  are  mutually  orthogonal  in  an  integral  sense  on  the  unit  sphere. 
The  spherical  scalar  wave  functions  are  defined  as 

where  fnl\'yv)  =  in(ir)  >  /n2^( 7r)  =  kn(^r)  are  modified  Bessel  functions. 
They  satisfy  the  Wronskian  relation 

W{si„(s),sA:n(s)}  =  stn(«)[sfcw(s)]'  -  [sr„(s)]'sfcn(s)  =  -1 

7  =  [ sfi(a  +  se)]1,  2  with  a,  e  are  permeability,  conductivity,  permittivity, 
respectively,  s  is  the  variable  of  the  two-sided  Laplace  transformation.  Co¬ 
efficients  times  the  scalar  wave  function  En,m,p(7C)  when  summed  over  all 
possible  indices  satisfy  the  scalar  wave  equation  which  for  each  function  we 
can  write  in  operator  form  as 

[V2  -  72]H!I.’m,(7>rl  =  0 
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From  the  solution  of  the  scalar  wave  equation  one  constructs  as  usual  the 
solutions  of  the  vector  wave  equation  of  three  kinds. 


1^)  — 

Mnln,p( 70  =  ^  X  [r^n,7n,p( 70] 

Nnjn,p( 70  =  X  ^n,m,p(7F) 

Note  that  all  three  kinds  of  vector  wave  functions  satisfy  the  vector  wave 
equation  in  Laplacian  form  which  we  can  summarize  as 


[V2  -  72 


f(D 

Lin^m^p 

M$n,p 

N{n}rn,p 


=  0 


We  can  also  write  a  curl  curl  wave  equation  for  only  the  second  and  third 
kinds  of  vector  wave  functions  as 


ivxv+^{ft}=° 

The  three  kinds  of  vector  wave  functions  have  some  interrelations  as 

^n,m,p( 7*0  =  -7T  X  tiPmAlf) 

=  -±V  X  N^m^r) 

A'','njn,p (70  ~  Mnjn,p{lfT) 

It  is  also  useful  to  write  them  as 

iJSnvM  =  [f^MYPn^  *)  +  [SinM]Qn,mA^  +)/V 
N$m,pM  =  +  hrIn)(V')]'Qn,m,p((,^)}hr 
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Plane  wave  in  spherical  coordinates 

As  shown  in  figure  3,  the  delta  function  plane  waves  (transformed)  can  be 
written  as 


f2e-7l>- 

f3e-7hr 


Y1  Z  Z  7^  +  &n,m,p 

n=l  m=0  p=e,o 

OO  71 

EEE  [^71,771,p'^71,771,p(/TO  ^7t,77l,jy 

7i=l  771=0  p=e,o 


n£L#M\ 


where 

'  ro  _  1  1/  iUt+1  2 7i +1  (n-77i)!  P^m)(cos(fli))  f  “  sin(m0i)  1 

njn,p  [  0,771  J v  J  7i(n+l)  (n+m)!"  sin(^i)  |  COS(m0i)  J 

t/  __  To  _  1  1/  i  \ 7i  2ti+1  (n-rn)\  dPjT\c os(6>i))  f  008(7/10!)  1 

n,m,p  L  0,77i J V  )  7i(ti-K)  (71+771)!  dOi  ^  sill(7M0i)  J 

Note  that  we  have 

-V  x  [r2e^f,'r1  =  f3e"7l"i  r" 

7 

-V  x  [f3e-7fi  r1  =  -l2e-^v? 

7 


which  is  associated  with  the  curl  relations  between  the  Mnln,p  and  Nn}n,p 
functions.  Furthermore  any  divergenceless  electric  field  expansion  (E)  can 
be  converted  to  a  magnetic  field  expansion  (//)  by  dividing  by  the  wave 
impedance  Z  of  the  medium  and  changing  Mnjn,p  to  -Nn}m,p  and  Nn\n,v  to 
Klnjn,p-  To  go  from  ll  to  E  multiply  by  Z  and  change  M,[ln4,  to  1 and 
yy(J)  to  —M^ 

yn,rri,p  lvl  n,rn,P’ 


Solution  of  the  scattered  field 

Define  our  incident  plane  wave  as  an  E  wave  (TM  wave) 


Einc(r,s)  =  E0l2e  7h'r 

f'hnc(r,s)  =  |^lV7flT 
^0 
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Expand  the  fields  for  r  <  a  as 


Ein(r,  s)  =  E{)J2  EE  KrnJCXA  V)  +  Cm,P^!U^] 

n=l  m=0  p=e,o 
rj  oo  n 

L(f,  s)  =  ^  £  E  £  iCjtfSwW  -  <•» 

^  n=lm=0p=e,o 

The  solution  of  the  scattered  fields  for  r  >  a  can  be  written  as 

°0  Tl 

L(f,  s)  =  ^oEEE  [<m,P^,P(70  +  C,m,pAtl,P(70] 

n=l  m— 0  p=e,o 
rp  co  n 

/?,c(r,  ^  S  S  [6n,m,PA/nfm,P(70  ~  <m,p^lp(7r1] 

n=l  m=0  p—e,o 

3.3  Perfectly  conducting  sphere 

Carl  Bamn  showed  in  [2]  that  there  only  exist  simple  poles  for  perfectly 
conductor  sphere.  Here  we  will  just  briefly  repeat  the  same  argument  in  our 
context.  Constrain  the  tangential  electric  field  to  be  zero  on  r  =  a,  we  have 
1 r  x  [Einc(r,  s )  +  Esc(f ,  s)]  =  0.  The  we  get 

C  X  [rt'„,m,P^,(.!m,P(70lr)  +  <m,pA4^.p(7«lr)]  =  0 
t  x  [6U,P^,p(7«t)  +  ^.P^,P(7alr)j  =  0 

This  give  equations  for  the  coefficient  as 


in(l/a)  , 


n’m’p  Jkn(7a)  "nm'p 


njri,p 


[7fl^n(7fl)]/  u 
[7afcn(7o)]'  n'm’p 


To  see  that  the  poles  must  be  simple  poles  we  just  need  to  show  that  all 
the  zeros  of  kn(s)  and  skn(s)  are  simple  zeros.  Since  kn(s)  satisfies  spherical 
Bessel  function  we  have 

s2-^kn(s)  +  2s-^j-kn(s)  -  [s2  +  7i (n  +  l)]fcn(s)  =  0 
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Suppose  the  zero  is  higher  than  first  order  say  a  2nd  order  zero  at  sa  ^  0. 
Since  both  kn(s)  and  k'n(s)  have  to  be  zero  at  sQ,  so  does  k'T[(s).  Thus,  the 
zero  must  be  at  least  a  third  order  zero.  Repeat  the  same  process,  we  will 
eventually  have  all  the  derivatives  at  sa  to  be  zero,  thus  the  function  must 
be  identically  zero.  So  there  exist  only  simple  poles  for  a"'m  .  The  argument 
for  bf”rnp  is  similar  as  [sfcn(s)]  satisfies  the  Riccati-Bessel  equation 


.s2  (P 
s2  +  n(n  +  1)  ds2 


Wi°Mi 


-  */?’(»)  =  0 


For  more  details  of  the  perfectly  conducting  sphere  including  surface  current 
and  charge  densities  please  see  [2]. 


3.4  Scatterer  of  spherical  sheet  impedance  loading 

3.4.1  Lossless  sheet  impedance 

Spherical  coordinates  (r,  0 , 0)  as  in  figure  3  are  one  of  the  few  coordinate 
systems  in  which  solutions  of  Maxwell’s  equations  are  separable.  In  partic¬ 
ular  let  us  assume  a  sheet  impedance  Zs(s)  (a  scalar)  which  is  located  on  a 
spherical  surface  give  by  r  —  a  and  which  is  independent  of  0 , 0.  This  sheet 
impedance  relates  tangential  electric  field  and  surface  current  density  as  in 
[3],  we  have 


17  •  E(a,  0 ,  <j>,  s)  =  Zg(s)Ja(6,  <j),  .s) 

It  —  1  —  lrh  =  transverse  dyad 
1  =  identity  dyad 

~  stands  for  two-sided  Laplace  transform.  The  surface  current  density  is  in 
turn  related  to  the  magnetic  field  via 

lr  x  [// («+,  0,  (j),  s)  -  II (a-,  0 ,  (j),  .s)]  =  0,  .s) 

The  sheet  impedance  function  Zs(s)  also  has  to  satisfy  Foster’s  Theorem  to 
guarantee  lossless  boundary  conditions. 
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Foster  Theorem 

In  1],  a  positive  real  function  F(s)  is  an  analytic  function  of  the  complex 
variable  s  =  a  +  ju,  which  has  the  following  properties: 

1. F(s)  is  regular  for  a  >  0 

2. F(a)  is  real 

3. a  >  0  implies  5R[F(s)]  >  0 

A  reactance  function  is  a  positive  real  function  that  maps  the  imaginary 
axis  into  the  imaginary  axis. 

Theorem:  A  real  rational  function  of  s  is  a  reactance  function  if  and  only 
if  all  of  its  poles  and  zeros  are  simple,  lie  on  the  ju- axis,  and  alternate  with 
each  other.  In  other  words 

X  =  ^2  +  ^)(6-2  +  ^)---(,2  +  u;22n_1) 

(s*  +  +  U1)  ■  ■  ■  (s*  +  ul) 

is  a  reactance  function,  where  k  =  2n  —  2  or  2 n,  K  >  0,  0  <  ljo  <  <^i  < 

*  *  *  <  ^2u-l  <  <^2n  <  OO 

Theorem:  A  rational  function  of  5  is  a  reactance  function  if  and  only  if 
it  is  the  driving-point  impedance  or  admittance  of  a  lossless  network. 

3.4.2  Solving  the  scattering  problem 

Matching  the  boundary  condition  on  r  =  a,  with  the  sheet  impedance  and 
continuity  of  the  tangential  electric  field  gives 

It  •  [Ei7ic(a+,d,<t>,s)  +  Eac(a+,d. </»,«)]  =  1,  •  Ein(a-,0,<p,s) 

=  Zs(s)js(0,<p,s) 

=  Zs(s)  x  ,  0,  (f)y  s)  “f-  H sc{^d~{- .  0)  (p,  s)  Hm(o  ,  6.  (f),  s)] 

Plugging  in  the  expansion  we  derive  a  system  of  equations  involving  a"  rn  , 

<m,p>  S°lve  f°F  an,m.p  and  6"  m<p  We  get 

//  _ 

a"'m’P  “  1  +  ^(7«)2*n(7«)M7a) 

i //  _  _ ^n,ro,p _ 

"'m'P  ”  1  “  ^j[7ai„(7o)]l7«fc„(7«)]/ 
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Now  the  surface  current  density  is 


Js(6,(f>,s)  =  -jr—  lt  •  R(a,  6,  <j>,  s) 
Zs{s) 


7  /  \  /  ,  /  ,  /  4  lan,m1pln('ya)  Rn,m,p(Q,  4>)  +  bn  m  p  Qn 

5W/  n=l  m=0  p=e,o 


3.4.3  Existence  of  2nd  order  poles 

The  coefficients  we  care  about  concerning  the  existence  of  a  second  pole  are 
Cl  =  l^an ,m,p  and  c2  =  Let’s  §ive  a  simple  example  showing 

that  a  second  pole  does  exist  for  coefficient  c2.  For  simplicity,  let’s  assume 
ja  =  s.  Consider  the  sheet  impedance  function 


Zs(s)  = 


(\e2  +  e+  ^)s 

+  ^  +  i 


Clearly  Zs(s)  satisfies  Foster’s  Theorem  with  K  >  0  and  u0  >  0.  The 
expansions  of  in(s)  and  kn(s)  are 

,  ,  ,  _  C  I  J)!2  's  ' 

"('Sj  «  j\  (n  —  j)\ 

in{s)  =  ^[(-l)n+1A:„(s)  -kn(-s)] 

For  n  =  0,  the  denominator  of  c2  is 

De{s)  =  (4 e~2s  +  4)  s2  +  (4e2  +  8e  +  4)s  +  1  +  2 e1-2*  +  2 e  +  e'2* 

It  is  easy  to  see  that  De(  —  ~)  =  0  and  ^De(s) |a=_i  =  0  or  in  Taylor  expan¬ 
sion  around  —  \ 

De(s)  =  ((16c  +  4  +  4e2)  (s  +  ±)V-f  e  -  |  e2)  (s  +  ±V+0  (  (s  + 


Thus  we  derive  a  second  order  pole  at  s  =  — 


0)] 
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In  general,  we  want  to  construct  a  sheet  impedance  function  Zs(s)  = 
such  that  c2  have  a  second  order  pole  in  the  left  half  plane  of  s  mean 
while  K  >  0  and  u  >  0.  The  denominator  of  r2  has  the  following  form 

De(s)  =  Ks+  s2in  (s) 

+  s^in  (s) 

We  want  to  solve  De(s)  =  0  and  -j^De(s)  —  0  for  /\,  uj  in  terms  of  s.  The 
solution  sQ  must  satisfy  sa  <  0,  f\(sa)  >  0  and  o;(sQ)  >  0.  For  n  =  0,  that 
is  to  solve 


-  s3-j-i,t  (s)  -  uj  in  (s)  -  u  s~~in  (s)  ]  kn  (s) 

as  as  ) 

-  s*£in  (s)  -  SLjin  (s)  -  S2U  -fin  (s))  ^~ku  (s) 

as  as  /  ds 


2  I\s  +  s2e  2s  +  s2  +  ue  2s  +  w  =  0 
2  K  +  2  .sr"2s  -  2  .sV2s  +  2  .s  -  2  w  c -2a  -  0 


The  solutions  are 

}<  =  _  s  (g~4^  +  2  e~2f)  T  1) 

2  se-25  4-  e~2s  +  1 
s 2  (— e~2s  +  2  se~2s  -  1) 
2  se~2s  +  e~2s  +  1 


From  figure  5,0.  approximately  when  s  is  chosen  from  —0.64  to  0,  both  K 
and  lj  will  be  positive.  Pushing  the  poles  to  even  a  higher  order  is  not  done 
here.  In  order  to  construct  a  high  order  pole  (including  the  2nd  order  case),  a 
transcendental  equation  lias  to  be  solved  analytically  which  in  general  is  not 
possible.  This  is  different  from  the  perfectly  conductor  sphere  case,  where 
only  a  system  of  linear  equations  need  to  be  solved. 

It  seems  to  us  that  it  only  works  for  coefficient  c2  with  n  =  0.  For  n  >  0, 
there  is  no  region  in  the  left  half  plane  of  s.  Although  we  don’t  have  a  rigorous 
proof  of  that,  it  seems  to  be  the  case.  We  tested  with  many  different  n,  A 
arid  uj.  K  and  uj  will  have  either  different  signs  or  both  will  be  negative. 
For  coefficient  cj,  it  doesn't  work  either.  We  tried  to  include  more  terms  in 
the  Zs(s)  according  to  Foster’s  Theorem  (i.e.  more  uj{)  with  larger  n,  but  it 
is  not  helpful  for  this  case.  Figure  7,8  show  some  results  of  different  cases 
with  different  n  ,  K  and  uji.  There  are  no  regions  for  which  I\  and  lj  are 
positive  simultaneously.  It  is  likely  that  for  lossless  sheet  impedance  loading 
boundary  condition,  most  scattering  poles  are  first  order. 
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Figure  5 


Figure  6: 
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Figure  7: 
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Figure  8: 


to  (cl  n=2) 
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4  Other  scattering  results 

Iii  [6],  a  canonical  problem  consisting  of  a  junction  in  a  transmission-line 
(single  conductor  plus  reference)  of  a  common  characteristic  impedance  Zv 
on  both  sides  of  the  junction  has  been  considered.  Scattering  from  such  a 
junction  is  only  forward  and  backward.  It  is  shown  that  scattering  from  a 
lossless,  reciprocal  network  at  a  transmission-line  junction  can  have  a  sec¬ 
ond  order  s-plane  pole.  In  [5],  Carl  extend  the  previous  results  to  a  more 
complicated  scattering  at  a  junction  of  transmission  lines.  A  transmission¬ 
line  structure  which  mimics  polarization  has  been  considered  to  make  the 
problem  even  closer  to  the  electromagnetic  case.  An  incident  wave  is  propa¬ 
gated  as  two  orthogonal  polarizations  on  a  four  wire  transmission  line.  This, 
in  turn,  scatters  forward  into  a  similar  four-wire  transmission  line  support¬ 
ing  two  orthogonal  polarizations.  Those  work  gives  an  closer  simulation  of 
the  electromagnetic-scattering  case,  contribute  to  the  discussion  concerning 
3-dimensional  electromagnetic  scattering  from  lossless,  as  well  as  perfect!} 
conducting  targets. 


5  Some  analysis  on  the  forward-scattering  the¬ 
orem  on  the  scalar  wave  equation 

5.1  Formulation  of  the  scattering  problem  from  a  scalar 
wave  equation 

We  start  from  the  simple  scalar  wave  equation  in  a  3-D  region  fi,  and  suppose 
c  =  1 


utt  =  V2u 

a(x)ut  +  j3(x)un  +  7  (x)u  =  0 

Let’s  call  Bu  =  0  the  boundary  condition,  un  is  the  outward  normal. 
Simply  applying  a  Laplace  transform  we  get 

s2u  =  V2n  Bit  =  0 

From  now  on,  all  the  computations  will  be  done  in  the  Laplace  s  domain, 
except  when  otherwise  indicated.  We  want  to  explicitly  compute  the  scat¬ 
tered  solution  u^sc\r}  ,s),  assuming  an  incident  plane  wave  u^inc\r,  ,s)  coining 
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in  the  ^-direction  (u  =  u^sc)  +  n(mc)).  Let  S\  denote  the  unit  sphere,  pa¬ 
rameterized  by  the  usual  angles  9 ,  (/?.  Suppose  u  is  the  incident  direction  . 
u  €  Si 


jJi  =  sin  ff  cos  (/?' 

<  uj2  =  sin  0 '  sin  ^ 

^  jJz  —  cos  8' 

Suppose  the  incident  plane  wave  in  the  time-domain  is  u^nc\x,t)  = 
A(L  —  u  •  x)  propagating  in  the  u  direction.  Then  u(mc)(f}  s)  =  /t(s)  •  e_&>u;  x 
(  just  the  Laplace  transform  of  u(mc)  ).  The  boundary  condition  becomes 
B(s)u(3C\x,  s)  =  —  B(s)u^tnc\x}  s)  ,  B  is  the  operator  for  the  boundary  con¬ 
dition.  Suppose  we  write  the  scattered  solution  u(sc\x,  s)  =  /l(s)  ■  U (x,  s\  u). 
Then  B($)U(x,s;u)  =  —  B(s)e~su>'x  on  the  boundary.  We  can  expand  U  in 
spherical  harmonics  outside  of  a  sphere  containing 

oo  n 

0 (x,  s; u)  =  Umn(s;u)Y™(9><p)kn(rs) 

n=0rn=— « 


When  r  — ►  oo,  A’n(rs)  ~  | ,  kn  is  the  modified  spherical  Bessel’s  function 
(good  at  infinity),  let  x  =  rq  ,  r  >  0  ,  q  €  S\  q  is  parameterized  by  (0,  y?), 
thus  £7 (£,  s;  a;)  =  U (r,  0,  </?,  s;  9\  ipf)  ~  (0,  </?,  0'.  </?',  s)  and 


*>,  (T,  .S)  =  ]T  J2  S~lUrn „(*\  ,  8)Y?{6,  <p)  =  v(q. 

n=0m= — n 


;,.s) 


5.2  Representation  Theorem 


-(.•sc) 


W  -  "5 


-(sc) 


(y)( 


3 


-|x-y|s  -|x-y|.s  Q 

■)  -  ti — —  t « sc  (y))]rfA 


3^|x-j/|  |x  —  y\  dvy 


Here,  the  parameters  :r,  y  are  vectors,  g^-is  the  outward  normal  derivative, 
cj,  s  have  been  omitted  for  simplicity  (always  there),  dAy  means  integrating 
over  y-area. 

Proof: 


V2 


e 


—rs 


r 


d2  2  d  e~rs 

3r2  r  di  r 


r 
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Figure  9: 


a 


D 


R 


Ball  of  radius  R  centered  at  x 
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r  is  the  radius  centered  at  x 


/  /  u(sc)V2?— -  {V2u{sc))(—dv  =  {) 


bv  Green’s  identity,  the  left  hand  side  also  equals  to 

d  e”lx~yb 

) - V 

2/1 


r  r  ft  p-\x-y\s  p-\x-y\s  ft 

/  S[i  5  W))^'+ 


(22) 


/  /  (23) 

\x-y\=R 

r  r  p-ts  p-ts  p-cs  £ 

j  J  4K2[i,“'(!/)(-s—  -  — )  -  ~(g~i(’‘\y))\dA,  (24) 

|x-y|=e 

dAy  means  integrate  over  a  sphere  with  surface  area  =  1.  In  equation(22) 
^ir(^(5c)(y))  ~  svi80) £-j^-  +  o(-^r)  and  it’s  easy  to  show  that  (23)  — >  0 

as  R  — >  00  for  3 R$  >  0.  Equation  (24)  — ►  4nu^(x)  as  e  — >  0  (behaves  like  a 
5  function). 


Thus,  we  derived  the  representation  for 

d  e“|x“y|‘s 


n^\x)  =  -~J  J  [«(sc)(y)( 


e-k-vl*  Q  ,  v 
)-TT— T(^-«(SC)(2/))]^y 


00 


|x  -  y I  |x  -  y I  <9^y 


For  the  far  field  pattern,  we  want  to  look  at  terms  like  limrerau^ac\x)  , 

r — >00 

because  that’s  what  v  looks  like.  (Note  fixed  origin  x  =  rx,  x  €  <Si) 


y|2  =  b’i:  -  t/p 


2rx  •  y  +  \y\2 


2, 


x  -y\  =  r\  1 - x  ■  y  +  =  r  -  x  •  y  +  O(-) 


e-\x-y\s  e(x*y)s+o(^) 


re 


I  r  —  2/1  1+0(7) 


e^x  y*s  as  r  —>  oo 
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Thus  ,  we  have  the  representation  for  v  in  the  far  field 


an 

Note  e^x  y)s  =  u(,nc)(y;  — i) 

Apply  Green’s  identities  inside  (for  incident  wave,  good  at  the  origin) 


0  =  JI  /«,"“i(j/;-;iV2u,,na(.v;  -.t)  ■  -x)dv, 

-  -i)]dAy 
on 

similarly,  apply  outside  il  and  use  the  radiation  condition  (5R,s  >  0) 


//> 

on 


du(sc\y;—x)  du^sc\y;u) 


dvv 


)-(- 


dVy 


-)u(sc)(t/;  -x)}dAy  -  0 


5.3  Reciprocity  from  boundary  conditions 

Let  us  assume  general  boundary  conditions  like 


dii 

a(rc)su  +  fl(x)- - h  7(x)u  =  0 

(si'll 


rewrite  as  +  A (s)u^  =  —  -  —  A (s)u^tnc) 

Consider  2  cases: 

i)  Dirichlet:  /?  =  0,  A  =  1 

ii)  mixed:  0  =  1;  ?RA(s)  >  0 


case  i): 
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v(.T,Ld.  S ) 


case  ii) 
v(x,u,s) 


-  -h I 

di} 

=  ~lj  /[«(“>(K*)(- 

an 


dit{inc){y;  -x)  du^iy-Sj)  ,itu 


{i’u)(y;  -x)]d' 


-ill  -*)  -  ~Ad.\ 


di/„ 


Oil 


rOu^(y,  -x). 


Oii^inc^ (y;  x) . 


-  -ill - **“•<* 

an 

=  ~H /t4W(K-i) 


9u(inc)(y;  u))  9?/(flC)(T/;  — £)  ,-nr 


an 

a(— ca,  5) 


1  _d^"c>(y;-.i),  du<sc>(y;£)..(jnc)  , 


4»J  / [i'  <K“)<  ft, 


;)- 


d«/„ 


-V'nc)(y\  -i)\ 


dil 


an 


+A(5)u(i"c)(y,  -x)[u(inc)(y;^)  +  «(sc)(y;^)] 

— A(s)fi^sc*(y,a;)[f^,,ic^(j/;  — x)  +  u^sc\y,  — x)] 

,  du('tnc'1  (y ;  —  x)  (y;  u>) 


-iff  (*""%;"><- 


') 


“<s,)(y;  -*)] 


an 


+A(s)[u<i"c>(y,  -x)u(,nc)(y;  w)  -  «<ac) (y, d.)ti(*c) (y;  -x)] 


-  -w 


du(mc)(y;w)  ,(sc)  ,(inc)  du(sc)(y;  -x) , 


3^, 


-u(sc)(y;-x)  -  u(!nc)(y;^)(- 


an 

x,  s) 
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There  are  more  general  cases  in  references  [8]  [7]. 


5.4  Forward  scattering  theorem  based  on  energy  the¬ 
ory 

Define  the  total  energy  to  be 


Integrate  the  outer  flux  term  f  *  §7  over  time,  we  have 

5 


f 


du  du 
dt  dr 


(bv  Parseval  ,  assume  real  u) 

f”  /  /'  ,  3u(y,S)  r 

=LJ  Jsu(s)-*-=L 


su(y ,  s) 


9u{y, -s) 
dr 


su(y ,  s) ■ 


du{y,  -s) 

dr 


+ 


^  (25) 

'  j  (26) 

5 


+  1  I  su{sc)(y,  s) 
s 


(sc),.  du(^(y,-s) 


dr 


(27) 


+  /  I  ***(,..)  ■***<*— >  (28) 


dr 
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term  (28)  when  R  is  large  and  s  purely  imaginary, 


c  c^s 

s-R-v(y,u;s)){s—v(y,d>]  -s))dA 


- II 


v(y,LJ]s)v(y,u]  —s)dA 


Si 


term  (25)  =  0  ,  since  u<mc)  =  e~Rsu =>  R2s2J  Juj  •  y  =  0 

S\ 

term  (26)  +  (27)  ~ 


l  C  A 

R2  J  J  (se~ltsw  y  ■  '—^-v(y,  u>;  - s )  +  — -v( y,  u\  s)  ■  su>  ■  yeIUu)  y 
5, 

=  fls2 1  I e^1-^  ■  v(y,  u;  -s)  +  u  ■  w;  s) 

' 

Suppose  Cj  =  C3,  in  coordinates  c u  •  y  =  cos  0 


J  -vix^-s) 


2tt 


3 /?«(!— cos  #) 


t’($,  (/?;  —5)  sin  OdOdip 


IT  L'v(^-s)TecR*"",°)dv 

O? r 

0, -*)«“■- »(0, 0, -*)] 

Ptt  1 

—  [u(tt,  (),  -s)em"  -  u(0,  (),  — *)]  +  o(-) 
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// 


Si 


n7T 

e-«s(I-cosfl) 


n(0,  y?;  s)  cos  9  sin  9d9dy 


1  /*27T  ^74  j 

—  /  /  n(0,<z>;s)cos0— e-*»(i-«.®)d 

Jo  Jo  du 


2  7T  , 

/?5 


(),-s)e  +  ^(0.  (),  — .s)] 


1  /»27T  /»7T  J 

J  J  jQ[v(d,ws)cose]eRs(1-co*0)ddd<p 

07 r  1 

=  -^[w(^,(),s)e_2fls +w(0,();s)] +  o(-) 

While,  the  last  step  is  not  completely  obvious  here.  We've  used  the 
method  of  stationary  phase  to  determine  the  asymptotic  behavior  for  R  large. 

For  example, 


/  cos0)dO  ~  -^[u{9,  (p; -a)  cos  9} 

combine  (25)  ~  (28)  energy  flux  for  R  large  ,  we  get 

du(y,  - s ) 


2tt  I 

6=0 

[R\s\\ 

1/2 


,171-/4 


su{y, s) 


dr 


=  S2  J  J  v{y,u;s)v(y,Lj-- 

Si 


s)dAy 


+27rs[t>(— —s)e2Rs  -  v(u/,ur,  — s)] 
+2irs[v(—u},  lj;  s)e~2Rs  +  v(Cj,  lj;  ,s)] 

since  27r.$[w(— u,lj\  —s)t2Rs  +  v(—u},lj;s)c  2/?s]  is  odd  in  s,  so 
f  2i[s[v{—u,u’,  —s)e2Rs  +  v{—Cj,Cj',  s)e-2/?s]  =  0 

Energy  flux: 


s'  i  i  v(y, u>;  s)v(y,  u;  —s)dAy  +  2ns[v(u ,  u;  s )  —  v(u,  a);  —  s)] 

5! 

For  a  lossless  scatterers  the  net  energy  flux=0  i.e. 
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A  27T 

v{y,u;s)v(y,u;-s)dAy  =  — — [v(&,<2»;  s)  -  v(u>,u>; -s)]  (29) 


Si 


which  is  the  same  as  the  ones  in  [4]  analogous  to  (2.4)  (5.6)  (6.6). 

5.5  Numerical  results 

Suppose  we  have  the  following  wave  equation  and  Dirichlet  boundary  condi¬ 
tion. 


s2u  =  V2u 


it  =  0  .  r  =  1 

Assume  an  incoming  plane  wave  is  in  the  2-direction  i.e. 

£(iuc)  _  e-s(0,0  l)  (x.T/,2)  _  e~sz 

We  want  to  compute  the  exact  scattered  solution 


5.5.1  Background  information 

Legendre’s  polynomial  and  spherical  harmonics  use  Rodriguez’  for¬ 
mula  the  Legendre  polynomial  can  be  written  as 


Pn(x) 


Orthogonality  relation: 


1  dn 
2 nn\  dxn 


(x2  -  1) 


/_,  pAz)r,.Wx  =  ~K'„ 

If  the  problem  does  not  possess  azimuthal  symmetry,  we  will  need  asso¬ 
ciated  Legendre  functions.  For  m  =  0,  ±1,  ±2, . . . ,  ±71 
Associated  Leyendre  functions  are 


2’*n! 


dxn+m 
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Properties  and  Orthogonality  relation: 


(n  +  m) 


s: 


p,',i{x)p,r;i(x)dx  =  2 

n  v  '  ”  v  '  2n  +  1  (n  —  m)\  n 


Notes: 


1.  Pn(x)  forms  a  complete  orthogonal  set  for  the  expansion  of  functions 
of  the  variable  x 

2.  Pnm(cos(0))  •  elTn<p  forms  such  a  set  for  the  expansion  of  arbitrary  func¬ 
tions  on  the  surface  of  a  sphere. 

Spherical  harmonics  are  defined  as 


o*.  v)  =  J 2"  t 1 !" ,  ”.*|;c(c°sw)  ■  c”"*’ 


4n  ( n  +  m)! 


Orthogonality  relation: 


f  '  f  r*  — ? - 

/  Y,:n(e,<p)Y;?(9,v)dn=  /  Y™(e,<p)Yy  (e,<P)smed6d<P  =  6nn'6TI 

JAir  JO  Jo 

Modified  spherical  Bessel’s  functions  Hankers  symbol 

r(§  +  n  +  k) 


(n.  k)  = 


k\ra  +  n-k) 


Modified  spherical  Bessel's  functions  are  the  solution  to  the  equation 
zSuj  +  2^u;  —  \z}  +  ti(ti  ~f  1)]cj  =  0 


•»(*-)  = 

[e~ni"l2jn{ze1"l2)  argc  G  (— 7t.7t/2] 
ye3nm/2j^ze-3Tn/2^  arg  Z  G  (7r/2,  7 r] 
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Expansion  of  in{z) 


in(z)  =  (2z)  1[R(n+^,-z)ez -(-l)nR(n  +  ^z)e  z] 

Where 

R(n  +  \,z)  =  ]CL-=o(n  + 


knW 


^  n+  j  (2) 

+  i  fc)(2«)-fc 

k= 0 

7T  „  1 

R(n  +  2’Z) 


Notes: 

zn,  kn  are  the  particular  solutions  we  will  be  using  later  on.  in 
good  at  origin,  while  kn  behaves  good  at  infinity. 

Some  useful  expansions 


=  £(2n+l)[J-/„+j(;)]P„(cos0) 

n=0  ’ 

oo 

=  y^(2  n  +  l)in(z)Pn(cos9) 


71=0 


, , — z  cos  6 


OO  j - 

=  +  ^r^2:ln+^^(cOS()) 

n=.0  ' 

oo 

=  y~](2n  +  l)(-l)wzn(z)Pn(cosfl) 


1=0 


5.5.2  Computation  of  the  solution 

We  already  have  the  expansion  for  the  incident  wave 

oo 

u(nK>  =  g-.»  =  e-»rcos e  _  y^(2n+  l)(-l)nin(sr)Fn(cos6>) 

71=0 


behaves 
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The  general  scattered  solutioii(in  this  case)  has  the  form 


oc  m=n 

«(9C)  =  E  E  amn(8)kn(sr)YZ'(0,tp) 

n=0  m=— n 
oc 

-  ^an(s)kn(sr)Pn(cos  6) 

71  =  0 


The  second  line  comes  from  azimuthal  symmetry.  In  our  case,  the  system 
is  symmetric  about  thus  m  —  0  i.e.  no  m  dependence,  then  we  have  a 
simplified  form. 

Applying  the  boundary  condition,  at  r  =  1 


a(s)fi(*c>  + 13 


=  —a(s)u 


(inc) 


^inc) 


dn  dn 

n  is  the  outward  normal,  i.e.  r  in  our  case  (sphere).  We  get 


an(s)[a(s)k„(s)  +  s/3k'n(s)}  =  (2  n  +  l)(-l)"+1[a(s)i7,(s)  +  s/?i',(s)J 

/  x  =  (~l)w+1(2n  +  l)[a(>);K(s)  +  sfii^js)] 

U[  ’  a(s)kn(s)  +  sPk'n(s) 

Then 


~(sc) 


-E 


n=0 


(  — l)n+1(2n  +  l)[a(s)in(s)  +  s^,(s)] 
a(s)kn(s)  +  S/3k'n(s) 


kn(sr)Pn(cos6) 


Let’s  look  at  the  far  field  pattern  r  — >  oc,  kn(sr)  — ►  | (asymptotic 
behavior) 


- («c)  =  £  V  (~1)T"l'1(2rt  +  +  3#»(g)]  P  m 

30  r  2  “j  s[o(«)M«) +  «/?<(«)]  nl  ^ 

By  our  previous  definition 
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t;(x,c2>;  5) 


7T  (-  l)n+1(2w  +  l)[o-(5)i„(.s)  +  s/^t(s)] 
2  ho  s[a{s)kn{s)  +  sf3k'n(s )] 

-  ^  Un(s)Pn (cOS0) 


For  this  v,  (29)  should  hold. 

For  simplicity,  let’s  assume  a  Dirichlet  boundary  condition  i.e.  a(s)  = 
l,/3  =  0  to  compute  the  singularity  expansion  and  check  the  numerical  re¬ 
sults. 


5.5.3  Singularity  Expansions 

We  will  perform  a  singularity  expansion  as  we  did  in  [4],  finally  we  will  see 
it  holds  (numerically  verified).  By  using  the  expansion  for  in.  kn  we  have 


iVt(s)  = 


(— l)n+1(2n  +  l)fn(s) 


sA^n(s) 


-  (— l)n+1(2rc  +  1) 


(2s)  1  [R(n  +  —  s)es  —  (— l)nR(n  +  s)e  s] 

n(2s)~1e~s  R(n  +  i  s)s 


(2 n  +  1)  r(-l)n+1fl(n+  ~s)c2 s  +  1' 


7 r 


R(n  +  s)s 


SJ 


Consider 


(-l)"+ifl(n  +  |,-a)  _  Qjjs) 

R(n  +  l,s)s  Q$(s) 

It  is  easy  to  see  that  Q",  Q"  are  polynomials  of  degree  n,  n+l  respectively. 
From  the  property  of  iu .  kn,  we  can  conclude  (.)'{ .  do  not  share  zeros  and 
zeros  are  simple. 

Suppose  Q%(s’‘)  =  0  for  j  =  0, 1, . . .  n 

qm  a  q1 

<%(«) "  hs- s" 
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Cn  9M1 
1  iQt(Sj) 


Proof  of  Cj': 


Suppose  QT2l(s)  =  B  UU(s  -  «"),  we  multiply 
by  ( s  —  srjf)  and  set  s  =  s™,.  We  get 


C" 


—  \  n  _ i 

/-^j  —  0  s-s 


svn  _ 

K/j'  — 


Q\l{s?) 


QUfA 

BnU*Asf-si)  iQZW 

Thus,  we  have  the  singularity  expansion 


Vn(s)  = 


2  Tt  ~h  1 


7T 


n  y^n 

[E^ 

J= o 


c?  2s  l 

e2A  +  - 
s"  sJ 


=  ZVAs)Pn(cOSd ) 


(30) 


n=0 


5.5.4  Connection  with  the  forward-scattering  theorem 

In  [4].  Dr.  Carl  Banin  suggests  the  following  form  for  the  scattering  dyadic, 

A  (lr,  l$;s)  =  >  - cQ{-lr)cQ{h)  +  entire f unction 

s  sa 

nr 

(In  our  case,  it’s  just  (30))  and  derives  the  following  equation, 


which  we  are  going  to  verify  in  the  next  part.  In  order  to  derive  an 
analogue  of  the  above  equation  from  (29),  we  will  follow  the  same  procedure 
as  in  [4].  Since  there  is  an  infinite  sum  in  the  integral,  we  will  also  use  the 
orthogonality  to  simplify  the  equation. 

In  (29)  ,  let's  multiply  both  sides  by  s  —  s™,  and  take  the  limit  s  — >  s™, 
(s"  is  one  zero  of  Q%  ,  s "  ^  0)  ,  using  the  orthogonality  of  PTl  we  get 


7T 

2 


2  n  +  1 


cys"Pn{  cos 


7T 


e)Y,M~s2)Pn(cose) 

n= 0 


2? t  2n  +  1 

•sq  * 


n„2s? 


die 


P„  (cos  0) 
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Table  1:  Coefficient  Cj  for  n  —  1  5  j  =  0,  ...n 


ri 


-1 

-1 

-1 

-1 

-1 

2 

2.564e-016  -  3.4641i 

-5.1583  +  0.67542i 

2.2877  +  6.9041i 

8.3088  -  5.0428i 

2.564e-016  +  3.4641i 

-5.1583  -  0.67542i 

2.2877  -  6.9041i 

8.3088  +  5.0428i 

12.317 

-2.2877  -  30.596i 

111.19 

-2.2877  +  30.596i 

-62.905  4  12.51  li 

-62.905-  12.51  li 

Table  2:  roots  of  Q^is)  for  n  =  1  5  j  =  1,  ...n  ($q  =  0) 


n 


-1 

-1.5  4  0.86603i 

-1.8389  4  1.7544i 

-2.1038  4  2.6574i 

-2.3247  4  3.57 Li 

-1.5  -  0.86603i 

-1.8389  -  1.7544i 

-2.1038  -  2.6574i 

-2.3247  -  3.571  i 

-2.3222 

-2.8962  4  0.86723i 

-3.6467 

-2.8962  -  0.86723i 

-3.352  4  1.7427i 

-3.352  -  1.7427i 

5  f  v(-sl)[Pn(cOs0)}2  =  -~Pn 

^  J  4ir 


(cos  0) 


7r  2n  +  1 
2~ 


7T 


«  /n»n 

j 

/  -j  _  <?n  _ 

J=0  a 


~2s" 


~S'„ 


r2i t 


[P„  (cos  (9)]2  sin  BdOd 


-r 


=  ~%Pn{C  OS0) 

C'TI 

/A 


Simplifying  we  get 


E 


j= 0 


P'lTl 

_ri _ e-2*s 


—  —  S; 


(32) 


Finally,  this  is  the  equation  analogous  to  (31),  which  we  are  going  to 
verify  numerically. 

CJ  are  the  coefficients  of  the  singularity  expansion,  sj  are  the  zeros  of 
Qj(s),  which  are  also  the  zeros  of  kn(s),  s™  is  anyone  of  s",  s7^  ^  0 

For  your  reference,  I  computed  all  the  CJ  ,  s",  error=|LHS-RHS|  of  (32) 
See  Table  1,  2,  3 
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Table  3:  Error  for  each  s™  and  for  n  =  1, 5 


n 


0 

1.1525e-014 

2.1245e-013 

4.1905e-012 

5.1391e-011 

1.2012e-014 

2.1245e-013 

4.1905e-012 

5.1391e-011 

3.6948e-013 

1.9759e-011 

6.9491e-010 

1.9723e-011 

3.9091e-010 

3.9096e-010 

5.5.5  Legendre’s  addition  theorem 

If  we  apply  Legendre’s  addition  theorem  [10]  to  the  above  results,  we  can  get 
a  more  general  solution  for  arbitrary  incidence  direction, 
let  7*i,  7*2  be  two  unit  vectors  with  components  given  by 

t* i  =  (sin  9\  cos  (j>\ ,  sin  9\  sin  (p\ ,  cos  9\ ) 

r2  =  (sin  92  cos  02,  sin  92  sin  02,  cos  92) 

The  angle  between  r\  and  f2  is  7  ,  defined  by 

r  1  •  f2  —  cos  7  =  sin  9\  sin  92  cos(<£i  —  (p2)  +  cos  9i  cos  92 


Pn{ cos 7)  =  Y"'  em— — ^|j-P™(cos  d\)P™ (cos  62)  cos  m(0i  -  <p2) 

(n  +  mV- 

to  =  1  and  em  =  2  for  m  >  0 

I11  terms  of  spherical  harmonics  ,  we  have 

'“(f‘  '  h)  =  2^Tl  £  OnWFTS) 

m=— n 

Thus,  the  incident  plane  wave  can  be  written  as 


'(znc) 


00  n  /  _  \  | 

V'(27i+l)(-l)njn(sr)  Y]  em 7— — r 7  PTn  (cos  #i )  P™ (cos  fl2)  cos  m 

(n  +  m)! 
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or  in  terms  of  spherical  harmonics 


in  which,  we  can  define  f\  as  the  incidence  direction  and  i'2  as  the  viewing 
direction  for  appropriate  choice  of  7. 

Thus,  the  scattered  solution(30)  becomes 


which  has  almost  the  same  form  as  appears  in  the  [4].  This  form  can 
be  used  to  determine  any  incident  direction  77  and  scattered  direction  for 
appropriate  choice  of  7. 

6  Conclusion 

We  show  that  for  acoustic  scattering  problems  arbitrary  order  scattering 
poles  can  be  constructed  with  certain  lossless  impedance  boundary  condition 
imposed  011  the  sphere.  While  for  the  hard  and  soft  sphere,  there  only  exist 
simple  poles.  For  the  electromagnetic  scattering  problems  for  spheres  with 
lossless  sheet  impedance  loading  boundary  conditions,  2nd  order  scattering 
poles  can  be  constructed.  There  only  exist  first  order  poles  for  the  perfectly 
conductor  sphere.  The  validity  of  the  acoustic  forward-scattering  theorem 
(scalar  wave  equation)  have  been  analyzed  as  well  as  some  numerical  exper¬ 
iment  of  the  theorem. 
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